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SYMAP1  -  AN  EXPERIMENTAL  SYMBOL  MANIPULATION  PROGRAM 


ABSTRACT 

SYMAP1  is  a  BRLESC  computer  program  designed  to  carry  out  a  variety 
of  algebraic  symbol  manipulations  including  arithmetic  operations, 
substitution,  certain  kinds  of  simplifications,  and  rudimentary 
differentiation.  Main  emphasis  has  been  on  polynomials  in  several 
variables  (including  truncated  power  series),  but  certain  other 
mathematical  forms  can  also  be  manipulated.  Examples  of  certain 
applications  to  numerical  analysis  and  theory  of  equations  are  included. 
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I.  INTRODUCTION  TO  SYMAP1 

In  order  to  relieve  a  researcher  of  certain  straightforward  but 
tedious  algebraic  symbol  manipulations ,  a  BRLESC  1  computer  program 
called  SYMAP1  has  been  developed  by  the  author.  Sufficient  capability 
has  been  demonstrated  to  permit  announcing  its  availability.  Additional 
features  will  be  added  as  they  are  found  useful,  and  more  powerful 
techniques  are  being  studied. 


Major  emphasis  in  the  early  stages  has  been 

on  manipulation  of 

polynomials  in  several  variables  (including  truncated  power  series), 

but 

many  features  of  SYMAP1  are  not  restricted  to  polynomials.  Principal 

manipulations  of  expressions  include  the  following  (see  also  Figure  l): 

Manipulations: 

Symbolic  Commands: 

1. 

Indicated  Siam  of  two  expressions 

ADD 

2. 

Indicated  difference 

SUBTR 

3. 

Indicated  product 

MULT 

4. 

Indicated  quotient 

DIV 

5- 

Indicated  power 

EXPON 

6. 

Certain  substitutions  in  expressions 

a.  Name  for  name 

SUBSTN 

b.  Name  for  expression 

SUBSTNE 

c.  Expression  for  name 

SUBSTE 

d.  Expression  for  expression 

SUBSTEE 

7- 

Display  of  any  previously  stored  labeled 

PRINT 

expression 


In  the  above  cases  the  expressions  have  few  restrictions.  In 
order  that  simplifications  of  the  results  may  be  virtually  automatic, 
however,  some  limitations  have  to  be  imposed  at  present.  Since  the  user 
may  want  to  retain  the  unsimplified  results  for  further  use, 
simplification  is  considered  a  separate  manipulation: 
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8.  Simplification  (limited) 

9.  Differentiation  (limited) 


SIMEXPR 

DIFF 


At  present  differentiation  is  limited  to  polynomial  -  like 
expressions  in  the  variable  of  differentiation  (negative  and  decimal 
fraction  exponents  are  allowed) .  Products  and  certain  more  general 
expressions  can  be  handled  by  considering  the  parts  separately  and 
combining  in  several  steps . 


Manipulations: 

Descriptions: 

ADD 

PA 

PB 

PCI 

Add  PA  and  PB  to  obtain  PCI 

SIMEXPR 

PCI 

PC 

and  simplify  as  PC 

SUBTR 

PA 

PB 

PD1 

Subtract  PB  from  PA  to  obtain 

SIMEXPR 

PD1 

PD 

PD1  and  simplify  as  PD 

MULT 

FT 

PA 

PEI 

Multiply  PT  times  PA  to  obtain 

SIMEXPR 

PEI 

PE 

PEI  and  simplify  as  PE 

DIV 

PB 

PT 

PF1 

Divide  PB  by  PT  to  obtain  PF1 

SIMEXPR 

PF1 

PF 

and  simplify  as  PF 

EXPON 

PT 

PK 

PG1 

Exponentiate  PT  to  PK  to 

SIMEXPR 

PG1 

PG 

obtain  PG1  and  simplify  as 

PG 

SUBSTEE 

PT 

X 

PA 

PHI 

Substitute  expression  PT  for 

SIMEXPR 

PHI 

PH 

expression  X  in  PA  to  obtain 
PHI  and  simplify  as  PH 

SUBSTEE 

PT 

Z 

QA 

QB1 

Substitute  expression  PT  for 

SIMEXPR 

QB1 

QB 

expression  Z  in  QA  to  obtain 
QB1  and  simplify  as  QB 

DIFF 

R1 

X 

DRl 

Differentiate  R1  w.r.t.  X  to 

obtain  DR1 


Figure  1.  Sample  Manipulations 
(Refer  also  to  Figure  2) 
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Given: 

PA 

1 . 3*X**2 .  - .  5*X**3  •  +6  .*Y*-*1 .  7 

PB 

X** -1 . +1 . 3*X**2 . - . 7*Y**1 . 7 

PI 

2.*Y**3. 

PK 

3- 

X 

X 

Y 

Y 

Z 

Z1+Z2 

QA 

(Zl+Z2)**2 . /Y+COS (XI ) +Z1+Z2 

R1 

3 .  *X**2 .*Y+X** -1 .  /  ( Y-tZ ) 

Results 

(which  can  be  manipulated  further): 

PCI 

PC 

l . 3*x**2 . - . 5*X**3 . +6 . *Y**l . 7+X** -l . +l . 3*X**2 . - . 7*y**1 . 7 

2 . 6*X+*2 .  - .  5*X**3  .  +X**  -1 .  +5 . 3*Y**1 . 7 

FD1 

PD 

1 . 3*X**2 . - . 5*X**3 • +6 . *Y**1 . 7- (X** -l . +1 . 3*X**2 . - . 7*Y**1 . 7 ) 
- . 5*X**3 . -X** -1 . +6 . 7*Y**1 . 7 

PEI 

PE 

2 .  *Y**3 *  ( 1  •  3*X**2 .  - .  5*X**3 .  +6  .■ *Y**1 . 7 ) 

2 . 6*X**2 . *Y**3 . -X**3 . *Y**3 .+12. *Y**H . 7 

PF1 

PF 

( X** -1 .  +1 . 3*X**2 .  - .  7*Y**1 . 7 )  /  ( 2 *Y**3  • ) 

.  5*x**-i .  *Y**-3 . + .  65*X**2.*Y*"*-3 .  - .  35*Y**  -l .  3 

PG1 

PG 

(2.*Y*-*3.)**3. 

8.*y**9. 

PHI 

PH 

l . 3* ( 2 .*Y**3 . )**2 . - . 5* ( 2 .*Y**3 . )**3 • +6 .*Y**l . 7 

6 .  *y**i . 7+5  •  2*y**6  .  -4 .  *Y**9  . 

031 

QB 

(2.*y**3.)**2./Y+C0S(XI)+2.*Y**3- 
2 .*Y**3 . +4 .*Y** 5 . +COS (XI ) 

DR1 

6 .  *X*Y  -X**  -2 .  *  (Y-tZ  )**-!. 

Figure  2.  Sample  Input  and  Results 
(Refer  also  to  Figure  l) 
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Certain  difficulties  sire  circumvented  by  doing  some  of  the 
manipulate  oris  on  algebraic  forms  expressed  in  a  non-standard  notation 
(called  here  canonical  notation) .  Provision  is  made  to  translate 
(simple)  expressions  to  and  from  this  canonical  notation  and  to 
manipulate  expressions  which  are  in  canonical  form: 


10. 

Convert  to  canonical  form 

CONVERT 

11. 

Reconvert  from  canonical  form 

RECONVERT 

12. 

Add  in  canonical  form 

CADD 

13. 

Subtract 

CSUBT 

ih. 

Multiply 

CMULT 

15- 

Divide 

CDIV 

16. 

Exponentiate 

CEXP 

IT. 

Simplify 

CSIME 

H 

00 

Differentiate 

CDIFF 

19- 

Substitute  expression  for  expression 

CSBEE 

20. 

Print,  one  term  per  line 

SPREAD 

The  last  (SPREAD),  usable  in  either  canonical  or  regular  notation, 
is  sometimes  more  convenient  than  PRINT  when  expressions  involve  many 
terms  in  several  variables. 

A  realistic  practical  example  will  be  considered  next,  using 
canonical  forms  largely. 

II.  EXAMPLE  FROM  NUMERICAL  ANALYSIS 

In  designing  numerical  approximations  for  partial  derivatives  for 
use  in  solving  partial  differential  equations,  the  following  expression 
can  be  considered  in  parabolic  cases: 

G(  I- ,T|)  =  (Ae-I=  +  B  +  Ce^Oe^  +  (De-^  +  E  +  Fe^)e"^1"' 

Here  A,  B,  C,  D,  E,  F  are  undetermined  weight  factors  to  apply  to 
function  values  at  points  of  a  rectangular  grid  neighboring  to  the 
point  at  which  the  derivatives  are  evaluated.  ^  and  T|  are  related  to 
step-sizes  in  two  orthogonal  directions  (grid  distances). 
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It  is  often  desirable  to  relate  the  step-sizes  by  setting 

2 

T1  =  P5 

where  P  is  another  undetermined  constant. 

Useful  relations  between  the  weights  A,...,F  can  be  obtained  by 

expanding  G(^,T|)  in  powers  of  |  and  T]  and  then  replacing  T|  or  by 

2 

expanding  instead  H(§)  =  G(5,P£  )  in  powers  of  §,  and  setting  the 
coefficient  of  each  low  power  of  5  equal  to  zero. 

The  currently  available  program.  SYMAP1  can  assist  in  this  effort. 

If  the  exponentials  are  replaced  by  truncated  power  series,  G(§,Tl) 
becomes  an  expression  involving  polynomials  in  5  an^-  T],  and-  H(§)  likewise 
in  ^  alone.  The  multiplying,  adding,  substituting  (in  some  cases), 
rearranging,  and  combining  of  like  terms  can  be  done  by  SYMAP1  on  the 
computer  BREESC  1. 

Figure  3  shows  a  possible  sequence  of  SfMAPl  steps  to  achieve  the 
result.  Truncation  of  all  series  through  a  given  power  of  §  is 
arranged  by  specifying  a  parameter  and  setting  a  program  option.  Terms 
through  power  8  can  be  handled  in  currently  available  memory.  Note 
that  G(5,T))  would  involve  several  dozen  terms  in  this  case,  a  situation 
in  which  hand  computation  might  easily  involve  careless  errors.  It 
was  to  avoid  such  errors  that  SYMAP1  was  originally  proposed. 

In  Figure  3  the  initial  input  (from  cards)  is  a  sequence  of 
labeled  expressions  which  are  to  be  operated  on.  The  end  of  this  set 
of  lines  is  indicated  by  a  special  sentinel  card  with  commas  in  card 
columns  1  and  2.  The  labels  are  arbitrary  (up  to  nine  characters)  but 
here  are  shown  as  rough  mnemonics  in  some  cases,  such  as  EXI  for  e^, 
EMHETA  for  e-^,  EHETX  for  e^  in  terms  of  etc.  The  power  series  are 

4  2 

truncated  here  after  |  or,  equivalently,  T)  . 

Following  the  sentinel  card  are  suitable  SYMAP1  instructions,  one 
per  card.  Many  alternative  sequences  could  have  been  used,  but  here  the 
labeled  expressions  wanted  are  all  converted  into  a  convenient  internal 
notation  called  canonical  notation  prior  to  any  combining. 
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Input : 


EXI  l. +XI**1 . +0 . 5*xi**2 . +. l666666?*XI**3 .  +  .  04l66667*XI**4 . 

emxi  l . -xi**l . +o . 5*XI**2 . - . i666666t*XI**3 . + . o4i66667*Xl**4 . 

EHETA  l.+0.5*ETA**l.+0.125*ETA**2. 

EMHETA  l.-0.5*ETA**l.-K).125*ETA**2. 


A 

B 

C 

D 

E 

r 


A**l. 

B**l. 

C**l. 

D**l. 

E**l. 

F**l. 


EHETX  1 . 40 . 5*P**  1 .  *XI**2 .  -K) .  125*P**  2 . *X  1**4 . 

EMHETX  1 . -0 . 5*P**1 . *XI**2 . +0 . 125*P**2 . *XI**4 . 

y  y 

Steps: 


CONVERT 

A 

CA 

CONVERT 

B 

CB 

CONVERT 

C 

CC 

CONVERT 

D 

CD 

CONVERT 

E 

CE 

CONVERT 

F 

CF 

CONVERT 

EXI 

CEXI 

CONVERT 

EMXI 

CEMXI 

CONVERT 

EHETA 

CEHETA 

CONVERT 

EMHETA 

CEMHETA 

CONVERT 

EHETX 

CEH1 

CONVERT 

EMHETX 

CEMHI 

CMULT 

CA 

CEMXI 

CMULT 

CC 

CEXI 

CADD 

T1 

CB 

CADD 

T3 

T2 

CMULT 

T4 

CEHETA 

CMULT 

CD 

CEMXI 

CMULT 

CF 

CEXI 

CADD 

Til 

CE 

CADD 

T13 

T12 

CMULT 

Ti4 

CEMHETA 

CADD 

SI 

S2 

RECO’IVERT 

S3 

GFCN 

SPREAD 

GFCN 

CMULT 

T4 

CEH1 

CMULT 

Tl4 

CEMHI 

CADD 

Sll 

S12 

RECONVERT 

S13 

HFCN 

SPREAD 

HFCN 

e« 


(*) 

-e 

e  w 

(*) 

T1 

e^2 

e4P?2 

A.e  ^ 

T2 

C.e^ 

T3 

Ae'^  +  B 

T4 

Ae~^  +  B  +  Ce^ 

SI 

(*) 

(Ae  ^  +  B  +  Ce^)e^ 

Til 

D.e‘5 

T12 

F.e^ 

T13 

De'?  +  E 

T14 

De  ^  +  E  +  F.e^ 

S2 

(*) 

(De"?  +  E  +  F .  e  ^ )  e 

S3 

(*) 

0(1,  Tl) 

(*) 

G(|,T1) 

(*) 

G(S,T1) 

511 

512 

513 

h(5)  =  g(5,p?2) 

(*)  only  if  G  itself  is  wanted 


Figure  3-  Input  and  Manipulations  for  First  Example 
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A 

+B 

+C 

+D 

+E 

+F 

-XI*A 

+XI*C 

-XI*D 

+XI*F 

+.5*XI**2.*P*A 
+.  5*XI**2.*P*B 
+.5*XI**2.*P*C 
- . 5*XI**2 .*P*D 
- . 5*XI**  2 . *P*E 
5*XI**2.*P*F 
+.5*XI**2.*A 
+. 5*XI**2.*C 
+.  5*XI**2.*D 
+.5*XI**2.*F 
- . 5*XI**3 • *P*A 
+.5*XI**3.*P*C 
+.5*XI**3-*P*D 
_.5*xi**3.*p*-f 
- . i66666t*xi**3 .*a 

+.i666667*xi**3-*c 
- . i666667*XI**3 .*d 
+ .  166666  7*XI**3  .*F 
+ .  25*XI**4 . *P*A 
+ .  25*XI**4  ,*P*C 
-.25*XI**4.*P*D 
-  .25*XI**4.*P*F 
+.125*XI**4.*P**2.*A 
+ .  125*XI**4  .*P**2 .*B 
+.125*XI**4.*P**2.*C 
+  .125*XI**4.*P**2.*D 
+  .125*XI**4.*P*-*2.*E 
+  .125*XI**4.*P**2.*F 
+.04l6667*Xl**4.*A 
+.04i6667*XI**4.*c 
+.04i6667*XI**4.*d 
+.04i6667*XI**4.*f 


Figure  4.  H(§)  in  "Spread"  Format 
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The  comments  at  the  right  indicate  the  result  obtained  at  each  step. 
(CADD,  CSUBT ,  and  CMULT  include  some  automatic  simplification.)  When  a 
result  of  sufficient  interest  is  obtained,  such  as  H(^)  near  the  end,  it 
can  be  reconverted  into  a  more  readable  form  and  if  desired  printed  one 
term  per  line,  i.e.,  "spread”. 

In  this  example  the  final  result  H(§)  is  displayed  in  ascending 
powers  of  Terms  involving  the  same  power  of  §  are  thus  adjacent  for 
the  convenience  of  the  reader.  (See  Figure  4.) 

Well-known  conclusions  that  can  be  drawn  from  this  display  are  that 
for  increasingly  higher  accuracy  of  numerical  approximation,  the  parameters 
A, . . . ,F,P  should  be  selected  to  satisfy  more  and  more  (from  the  top)  of 
the  following  relations,  one  for  each  power  of 


A+B+C+D+E+F  =0 

-A+C-D+F  =0 

P*(A  +B+C-D-E-F)  +(A+C+D+F)  =0 

3*P*(.  A+C+D-F)+(-  A+C-D+F)  =0 


6*P*(A  +  C  -  D  -  F)  +  3*P2*(A  +B+C+D+E+F)  +  (A  +  C  +  D  +  F)  =0 

Note  that  truncating  after  the  fourth  power  of  |  still  leads  to  42 
terms  in  the  expression.  This  number  would  be  much  greater  for  6th, 

8th,  etc.  powers,  with  increasing  chance  of  human  error  if  done  by  hand. 
Expansions  through  power  8  have  been  run  on  BRLESC  1  for  this  problem. 

In  Figure  3  actual  numerical  coefficients  (such  as  .125  and 
. 16666667 )  were  shown.  After  sufficient  manipulation  and  combination 
such  coefficients  might  become  unrecognizable;  so  they  could  be  treated 
as  additional  parameters  and  replaced  at  a  later  stage  if  desired.  (See 
Figure  5*)  This  would  make  each  term  longer,  however,  and  thus  require 
some  additional  machine  effort.  (There  is  a  limit  on  the  number  of 
factors  per  SYMAP1  term,  currently  seven  in  most  cases,  fourteen  in 
special  cases.)  Figure  5  also  further  demonstrates  that  negative 
exponents  of  simple  variables  can  be  handled  in  SYMAP1.  Decimal 
fractions  in  exponents  are  also  allowed. 


Given: 


EXI1  1 .  +XI**1 . -tTWO**  -1 .  *XI**  2 .  -HDWO**  -1 .  *TRE*-*  -1 .  *XI**3 . 

EXI2  TWO**-3.*TRE**-l.*XI**4. 

EMXI1  1 . -XI**1 . +TWO** -1 . *XI**  2 . -TWO** -1 . *TRE**-1 .*XI**3 • 

EMXI2  TWO**-3.*TRE**-l.*XI**4. 


EHETX1 

l.+TWO** 

-1 .  *P**  1 .  *XI**  2 .  -HTWO**  -3  .' *P**2  .' *XI**4 

EMHETX1 

1 .  -TWO** 

-1 .  *P**1 .  *XI**2 .  HdWO**  -3 .  *P**2 .  *XI**4 

A 

A**l. 

B 

B**l. 

C 

C**l. 

D 

D**l. 

E 

E**l. 

F 

F**l. 

TWO 

962 

TRE 

963 

CTWO 

2.$ 

CTRE 

>  > 

34 

Manipulations: 

Results 

CONVERT 

A 

CA 

CONVERT 

B 

CB 

CONVERT 

C 

CC 

CONVERT 

D 

CD 

CONVERT 

E 

CE 

CONVERT 

F 

CF 

CONVERT 

EXI1 

CXI 

CONVERT 

EXI2 

CX2 

CONVERT 

EMXI1 

CMX1 

CONVERT 

EMXI2 

CMX2 

CONVERT 

EHETX1 

CHI  e*P52 

CONVERT 

EMHETX1 

-^P|2 

CMH1  e  s  s 

Figure  5(a).  Modified  Version  of  First  Example.  Part  1 
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Manipulations : 


CADD 

CXI 

CX2 

CADD 

CMX1 

CMX2 

CMULT 

CA 

CEMX 

CMULT 

CC 

CEX 

CADD 

T21 

CB 

CADD 

T23 

T22 

CMULT 

CD 

CEMX 

CMULT 

CF 

CEX 

CADD 

TJl 

CE 

CADD 

T35 

T32 

CMULT 

T24 

CHI 

CMULT 

T34 

CMH1 

CADD 

S21 

S22 

RECONVERT 

S23 

HXI 

SPREAD 

HXI 

CSBEE 

CTWO 

TWO 

CSIME 

S31 

S32 

CSBEE 

CTRE 

TRE 

CSIME 

S34 

S3  5 

RECONVERT 

S3  5 

HI 

SPREAD 

HI 

>  > 


Results: 


CEX 

e^ 

CEMX 

e’S 

T21 

Ae_l= 

T22 

Ce^ 

T23 

Ae"^  +  B 

T24 

Ae-^  +  B  +  Ce^ 

T31 

De~  ~ 

T32 

Fe^ 

T33 

De“^  +  E 

T34 

De  ^  +  E  +  Fe^ 

S21 

S22 

S23 

H(S) 

H(5J  TWO, TRE) 

S23 

S31 

Replace  TWO 

S32 

S34 

Replace  TRE 

H(§),  simplified 

Figure  5(b).  Modified  Version  of  First  Example.  Part  2 
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III.  INTERNAL  NOTATION  FOR  CONSTANTS,  VARIABLES,  AND  FUNCTION  NAMES 

Variables  and  parameters  which  occur  in  the  manipulation  of  algebraic 
expressions  are  represented  in  the  preferred  internal  notation  of  SYMAP1 
by  three  digit  integers:  100,101, ... ,999.  Translations  to  and  from  this 
internal  notation  make  use  of  several  conversion  tables,  based  on  the 
length  of  the  external  symbol.  Within  the  computer  program  Table  1AV 
contains  single  character  names,  such  as  A,  B,  X,  Y  and  their  three 
digit  integer  synonyms.  Table  2AV  likewise  handles  two-character  names, 
such  as  XI,  NU,  Al,  X',  etc.  Tables  3AV,4AV, . . . ,SAV  are  provided  for 
names  of  three  or  more  characters. 

The  selection  of  the  three  digit  integers  is  arbitrary  but  does 
affect  the  lexicographic  sequencing  of  factors  within  terms  and  of  terms 
within  sums,  etc.  Thus,  if  some  final  result  is  wanted  in  ascending 
powers  of  Y,  say,  and  in  case  of  ties  in  ascending  powers  of  X,  then  the 
integer  for  Y  should  be  less  than  the  integer  for  X  and  both  should  be 
less  than  the  integers  of  other  variables  and  parameters. 

If  for  instance  we  represent  Y  by  501,  X  by  502,  and  A,  B,  C  by  601, 

2 

602,  603  respectively,  then  internally  Y  becomes  501**1.$  and  X  becomes 
501**2.$,  etc.  The  binary  representation  of  the  digits  and  characters 
within  BRIESC  1  and  the  ordering  rules  of  SYMAP1  then  cause  algebraic 
terms  in  Y,  X,  A,  B,  C  to  be  sequenced  as  follows: 

Y  ,YX  ,YXA,YXB,... 

y^,y^x,y2xa, . . . ,y^x^,y2x2a, . . . 

y5,y5x,y5xa, . . . ,y5x2,y5x2a, . . . 

•  •  • 

X  ,XA  ,XAB,XB, . . . 

2  2 

X  ,X  A,... 

X5,X5A, ... 

•  •  • 

Thus  terms  with  like  powers  of  Y  are  adjacent,  and  for  a  given  power  of 
Y  those  terms  with  like  powers  of  X  are  adjacent,  etc.  Numerical 
coefficients  do  not  affect  this  ordering.  (Sufficiently  high  and 
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negative  exponents  do  affect  the  ordering,  but  the  adjacency  within 
groups  is  maintained.)  This  adjacency  property  is  useful  in  combining 
"like"  terms  after  expansions. 

If  there  is  a  preferred  lexicographic  order,  this  should  be 
considered  and  the  entries  in  Tables  LAV,  etc,  specified  appropriately. 

Function  names,  similarly,  are  represented  by  integers  of  either  1 
or  2  digits,  using  the  same  Tables  IAV,  etc.  For  sequencing,  however, 
these  integers  are  considered  10000  times  larger  at  present,  in  order  to 
keep  functions  at  the  right  of  any  term  in  which  they  appear. 

Since  variables  and  function  names  are  represented  by  integers,  it 
is  necessary  that  actual  constants  be  distinguished  with  special  marks. 
In  SYMAP1  this  is  done  by  appending  a  dollar  sign  at  the  end  of  each 
constant.  Also  each  constant  is  required  to  contain  a  decimal  point. 
Thus  unity  is  l.$  in  canonical  notation  and  zero  is  .0$  internally. 
Decimal  fractions  at  present  are  carried  with  up  to  eight  digits  if 
positive  or  up  to  seven  digits  if  negative  (an  arbitrary  decision  but 
convenient  in  a  computer  with  10  characters  per  word) . 

Internally  each  factor  has  an  explicit  exponent,  unity  if  nothing 
else,  and  each  term  has  an  explicit  coefficient.  A  coefficient  or 
exponent  of  unity  is  removed  in  the  reconversion  process,  however,  to 
keep  most  displayed  expressions  as  near  to  normal  notation  as  possible. 

Let  us  consider  another  example  in  several  variables. 

IV.  EXAMPLE  FROM  THE  THEORY  OF  CUBIC  EQUATIONS 

As  another  example  of  symbol  manipulation  in  SYMAP1  let  us  consider 
the  nature  of  the  roots  of  a  cubic  equation.  Given  the  polynomial 

5  2 

p(x)  =  yi  +  bx  +  cx  +  d  , 

where  b,  c,  d  are  real,  let  us  characterize  the  region  of  bcd-space 
where  p(x)  has  three  real  roots.  It  is  well  known  that  the  requirement 
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for  this  condition  is  that  the  cubic  discriminant  defined  by 

A(b,c,d)  =  l8bcd  -  4b^d  +  b2c2  -  4c^  -  27d2 

be  non-negative.  Let  us  confirm  this  relation. 

Of  the  three  roots  of  any  cubic  with  real  coefficients,  at  least 
one  root  is  real  and  the  other  two  may  be  real  and  equal,  real  and 
unequal,  or  complex  conjugates.  Therefore,  we  can  let  the  roots  be  r, 
s+t,  and  s-t  with  r  and  s  real  and  t  zero,  positive,  or  a  pure 
imaginary  for  the  three  cases. 

The  real  coefficients  b,  c,  and  d  are  related  to  r,  s,  and  t  as 
follows : 

b  =  -  r  -  2s 

2  2 

c  =  2rs  +  s  -  t 
d  =  -  r(s2  -  t2) 

If  these  last  expressions  are  substituted  for  b,  c,  and  d  in  the 
discriminant  formula,  we  obtain: 

o  k  6 

A  =  Ao  +  Alt  +  A2t  +  A3t 

where  the  coefficients  A.  are  polynomials  in  r  and  s. 

If  t=0  p(x)  has  at  least  a  double  root  s  and  the  properties  of  the 

discriminant  state  that  A  vanishes.  Hence  we  note  that  it  must  follow 

that  the  coefficient  Aq  vanishes  identically: 

A  =  0  . 
o 

This  fact  can  be  verified  using  SYMAP1. 

If  indeed  A  =0,  then  our  problem  of  showing 
o 

A(r,s,t)  ^  0 

reduces  to  the  equivalent  problem  of  verifying  that 

;  2  2  4 

A/t  =  A^  +  A^t  +  A^t  , 

a  quadratic  in  t  ,  is  non-negative  (when  r,  s,  and  t  are  real) . 
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Input : 


B 

C 

D 

TM2 

TM4 

IT2 

n4 

IZERO 

14 

118 

127 


Steps: 

CONVERT 

CONVERT 

CONVERT 

CMULT 

CMULT 

CMULT 

CMULT 

CMULT 

CMULT 

CMULT 

CMULT 

CMULT 

CMULT 

CMULT 

CMULT 

CMULT 

CSUBT 

CADD 

CSUBT 

CSUBT 

RECONVERT 

SPREAD 


-l.*R-2.*S 
2.*R*S+S**2.-T**2. 
-1 . *R*S**2 . +R*T**2 . 
T**-2. 

T**-4. 

801**2.$ 

801**4.$ 

.0$ 

4.$ 

184 

274 


Results: 


B 

IB 

b 

C 

IC 

c 

D 

ID 

d 

ID 

ID 

ID2 

d2 

IC 

IC 

IC2 

2 

c 

IB 

IB 

IB  2 

b2 

IC2 

IC 

IC3 

c3 

IB2 

IB 

IB3 

b3 

IB 

IC 

IBC 

be 

IBC 

ID 

IBCD 

bed 

118 

IBCD 

T1 

l8bcd  =  term  1 

IB5 

ID 

IB3D 

b3d 

14 

IBJD 

T2 

4b3d  =  term  2 

IB2 

IC2 

T3 

2  2 

be  =  term  3 

14 

IC3 

t4 

4c  =  term  4 

127 

ID2 

T5 

27d2  =  term  5 

T1 

T2 

S2 

S2 

T3 

S3 

S5 

T4 

s4 

S4 

T5 

S5 

A 

S5 

DISC 

(No  term  vitho- 

DISC 

so  A  =  0) 

Figure  6(a).  Generating  the  Cubic  Discriminant 
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SYMAP1  can  be  used  to  determine  the  literal  form  of  the  coefficients 
A2,  and  A^  and  then  the  form  of  the  quadratic  discriminant 

E  =  (A2)2  -  4AjA?  . 


For  all  real  r  and  s  we  would  expect  to  find 

A^  ^  0  and  E  ^  0  . 

Let  us  verify  these  relations  also  using  SYMAP1  as  an  aid. 

Figure  6  indicates  a  possible  set  of  SYMAP1  steps  in  support  of 
such  verifications.  Input  forms  include  b,  c,  d  in  terms  of  r,  s,  t, 
four  useful  powers  of  t  (two  of  them  in  internal  or  canonical  form) , 
and  four  constants  (also  in  canonical  notation  for  convenience).  The 
three  digit  integers  801,  802,  and  803  were  selected  to  represent  t,  s, 
and  r  respectively.  Any  results  are  grouped  by  like  powers  of  t. 

Manipulations  include  converting  b,  c,  and  d  to  canonical  form, 
building  up  A  step  by  step,  displaying  A,  isolating  the  coefficients 
Ap  A2,  and  A^,  and  forming  and  displaying  the  quadratic  discriminant  E. 
The  observer  can  draw  his  conclusions  after  some  additional  hand  work 
including  some  factoring  not  within  the  capabilities  of  STMAP1. 

Figure  7  shows  the  literal  forms  of  A,  A^,  Ag,  A^,  and  E  in  "spread" 
form.  Note  in  particular  that  in  A  no  term  is  independent  of  t;  so 
Aq  =  0  as  expected.  A^  is  found  to  be  the  constant  4.  A2  can  be  factored 

by  hand  as  A.  =  -  8(s-r)2,  real  and  non-positive.  Similarly  A.  =  4(s-r)  , 

^  2  “*■ 
real  and  non-negative.  E  =  A2  -  4A^A^  =  0,  a  fact  which  might  not  be 

expected  without  further  analysis,  but  the  desired  conclusion  E  <  0  is 

certainly  verified. 

Replacing  A^,  A£,  A^  in  A/t  leads  to: 

A/t2  =  4(s-r)^  -  8(s-r)2t2  +  4t^ 

A  =  4t2[(s-r)^  -  2(s-r)2t2  +  t^] 

=  4t2[(s-r)2  -  t2]2  s  0 

We  have  thus  verified  that  the  discriminant  A  is  non-negative  if  the 
cubic  has  three  real  roots.  Furthermore, 
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CONVERT 

CONVERT 

CMULT 

TM2 

TM4 

ITM2 

ITM2 

ITML 

S5 

S5A 

t  c 

f4 

A.t-2 

RECONVERT 

S5A 

RS5A 

p 

I. 

SPREAD 

RS5A 

\  +  v  + 

V 

CSBEE 

IZERO 

it4 

S5A  A2A 

\  +  v2  + 

ZyO 

CSIME 

A2A 

A2 

Ai  +  A2t2 

CSUBT 

S5A 

A2 

A3 

v4 

CMULT 

ITML 

A3 

EEL3 

A  4-^ 

Ajt  .t 

=  A5 

CSBEE 

EERO 

it  2 

A2  AkA 

A1  +  A2.0 

CSIME 

AbA 

IDELl 

=  A1 

CSUBT 

A2 

EE  LI 

A5 

A2t2 

CMULT 

ITM2 

A5 

IDEL2 

A2t2.t’2 

<r 

11 

CMULT 

IDEL1 

EEL3 

Til 

V* 

CMULT 

IU 

Til 

T12 

CMULT 

IDEL2 

IDEL2 

T13 

A2 

CSUBT 

T13 

Til 

E 

A2  -  ^ 

=  E 

RECONVERT 

IDEL1 

DEE 

SPREAD 

DELI 

RECONVERT 

IDEL2 

DEL2 

SPREAD 

DEL2 

RECONVERT 

IDEL3 

DEL3 

SPREAD 

DEL3 

RECONVERT 

IE 

E 

SPREAD 

;  7 

E 

Figure  6(b) 

.  Analyzing  the  Cubic  Discriminant 
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A(R,S,T) 


-l6.*T**2.*S*R**3- 
+24 .*T**2 .*S**2 .*R**2 . 
-l6.*T**2.*S**3-*R 
+4.*T**2.*S**4. 

+4.*T*-*2.*R**4. 

+i6.*t**4.*S*r 
-8.*T#*4.*S**2. 

_8.*p**E.*R**2. 

+4.*t**6. 

-l6.*S*R**3-  MR,S) 

+24.*S**2.*R**2. 

-l6.*s**3.*R 

+4.*S**4. 

+4.*R**4. 

l6.*S*R  Ap(R,S) 

-8.*S**2. 

-8.*r**2. 

4.  a5(r,s) 


.0 


Figure  7.  Results  of  Second  Example 
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A  =  ht2[s-r-t]2[s-r+t]2 

=  l|-t2[(s-t)  -  r]2[(s+t)  -  r]2  s  0  . 

In  this  last  form  it  is  quite  clear  that  the  cubic  discriminant  A 
vanishes  if  t=0  or  if  r  =  s+t  or  if  r  =  s-t,  exactly  those  cases  when 
at  least  two  roots  are  equal,  and  A  does  not  vanish  in  any  other  case. 

In  summary, 

A  £  0  for  any  three  real  roots 
A  =  0  for  at  least  two  equal  roots  (all  real) 

A  >  0  for  three  distinct  real  roots 
and  A  <  0  in  other  cases,  namely  for  one  real  root  and  two 

complex  conjugate  roots. 

SfMAPl  could  be  used  in  a  similar  manner  to  help  derive  less 
well-known  relations  occurring  in  original  investigations. 

V.  SIMPLIFICATION  OF  FORMS 

Several  types  of  simplification  occur  automatically  when  SYMAP1  is 
used.  Others  can  be  achieved  by  appropriate  combinations  of  steps, 
specified  by  the  user.  More  flexibility  is  available  when  forms  are 
retained  in  canonical  notation;  so  use  of  this  notation  is  encouraged. 

In  general  the  following  simplifications  occur,  but  there  are 
exceptions,  especially  if  complex  combinations  of  steps  are  carried  out 
without  calling  for  intermediate  simplifications.  In  the  formulas 
below  a,  b,  c,  n,  m,  p,  q,  r  are  real  decimal  constants  (not  necessarily 
integers);  x,  y  are  primitives;  and  t,  u,  v,  w  are  terms  similar  to  the 
terms  of  a  polynomial. 

Within  a  term  like  factors  are  combined,  and  within  a  sum  like 
terms  are  combined.  Thus: 

n  m  p 

xx  -*  x  where  p  =  n  +  m 

nm,.nm  nm  ,  ,  , 

ax  y  +  bx  y  -*  cx  y  where  c  =  a  +  b 
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Powers  of  simple  products  are  expanded: 

(axnym)C  -*  bx^y^  where  b  =  aC,  p  =  c-n,  q  =  c>m 
The  distributive  law  is  applied  in  one  direction: 


axC(b-jXn  +  b2xm)  ->  r.jXP  +  r2xq 

but  the  reverse  process  (factoring)  is  not.  Special  cases  of  factoring 
can  be  done  in  several  steps,  however. 

Parentheses  are  removed  when  possible  in  special  cases,  such  as: 

t  +  (u+v-w)-*t+u  +  v-  w 
t  -  (u+v-w)-*t-u-v+w 

Simple  denominators  are  replaced: 

u/a  -*  b-u  where  b  =  a  ^ 

u/(axn)  -*  b.u-x  n  where  b  =  a  ^ 

u/(ax  +  by  )  -*  u-(ax  +  by  ) 


Special  cases  of  the  above  rules  include: 


n  -n 
x  *x 


n  n 
ax  -  ax  -*  0 


Indeterminate  cases  must  be  avoided  by  the  user: 

0°  0/0  0*  co  0°°  m°  °°/o°  1°°  00-00 

Resequencing  in  some  lexicographical  order  is  standard  to  group 
similar  factors  and  terms. 


VI.  DETAILS  OF  SPECIFICATION  OF  MANIPULATIONS 

The  principal  manipulation  types  were  introduced  earlier,  and 
examples  of  their  use  were  given.  For  those  interested  additional 
details  are  given  here . 

At  present  manipulations  must  be  specified  in  advance  on  formatted 
punched  cards.  (When  they  become  available,  on-line  keyboards,  light 
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pens,  etc.  will  make  SIMAP  even  more  responsive  to  user  needs.)  The 
manipulation  type  is  specified  at  card  column  1  and  any  additional 
parameters  at  columns  11,  21,  31,  ^1  as  needed.  The  parameters  are 
always  names  of  previously  specified  forms  or  variables  or  of  the  form 
currently  being  generated.  Up  to  9  letters  and  digits  may  be  used  in 
names . 

Arithmetic  operations  (add,  subtract,  multiply,  divide, 
exponentiate)  have  three  parameters:  two  operands  and  a  result,  in  the 
usual  order.  PRINT  and  SPREAD  require  only  one  parameter,  the  name  of 
the  form  to  be  displayed.  CONVERT,  RECONVERT,  and  the  simplify  types 
have  two  parameters  each:  the  given  form  and  a  (new)  name  for  the 
result.  (Names  cannot  be  reused  with  new  meanings,  even  if  the  old 
form  is  not  needed  further.) 

The  substitution  manipulations  have  four  parameters  in  the  order 
a,  b,  c,  d,  where  a  represents  the  form  substituted,  b  the  form  substi¬ 
tuted  for,  c  the  form  substituted  in,  and  d  the  result.  Thus 
CSBEE  CA  CB  CC  CD  means  "substitute  the  form  whose 
name  is  CA  for  the  form  whose  name  is  CB  wherever  it  occurs  in  the  form 
whose  name  is  CC  and  store  the  result  with  the  label  CD".  Type  CSBEE 
is  the  most  commonly  used  substitution  in  SYMAP1.  Type  SUBSTEE  is  very 
similar  but  uses  non-canonical  forms. 

SUBSTN  is  merely  a  renaming  operation:  "substitute  the  name  a 
(not  the  form  whose  name  is  a)  for  the  name  b  wherever  b  occurs  in  the 
form  called  c  and  label  the  result  as  d" .  SUBSTNE  (name  for  expression) 
is  similar  except  that  b  is  the  name  of  the  form  being  substituted  for. 
Finally  SUBSTE  (expression  for  name)  says  to  substitute  the  expression 
whose  name  is  a  for  the  name  b  wherever  b  occurs  in  c  and  call  the 
result  d.  All  of  these  operate  on  non-canonical  forms.  Simplification 
of  the  results  of  these  operations  with  SIMEXPR  is  not  always 
successful,  however;  so  use  of  canonical  forms  and  CSBEE  followed  by 
CSIME  is  recommended. 
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The  differentiation  operations  CDIFF  and  DIFF,  for  canonical  and 
non-canonical  forms,  respectively,  are  currently  limited  to  polynomial¬ 
like  expressions  in  the  variable  of  differentiation.  Extensions  are 
planned.  CDIFF  CA  Y  CB  means  "differentiate  the  form  whose 
name  is  CA  with  respect  to  the  primitive  Y  and  call  the  result  CB". 

Since  CA  is  canonical,  the  variable  Y  does  not  appear  in  CA  as  such  but 
rather  as  some  previously  assigned  three  digit  number  such  as  101. 

SYMAP1  looks  up  the  equivalent  101,  and  differentiates  with  respect  to 
that  symbol,  giving  a  canonical  result. 

Manipulations  are  carried  out  in  the  order  specified  by  the  sequence 
of  cards,  and  only  input  and  previously  completed  results  can  be  used 
in  later  manipulations .  The  final  operation  should  be  followed  by  a 
sentinel  card  (commas  in  columns  1  and  2)  to  terminate  manipulations. 

VII.  EVALUATION  OF  EXPRESSIONS 

Although  the  power  of  SYMAP1  lies  in  its  ability  to  manipulate 
polynomials  and  similar  non -numerical  expressions,  there  are  occasions 
when  an  algebraic  form  must  be  evaluated  for  particular  values  of  its 
variables.  This  capability  is  included  in  the  substitution  and 
simplification  operations  already  described. 

Given,  say,  a  polynomial  in  X  and  Y  where  the  values  X  =  4  and 
Y  =  -  5.25  are  to  be  assigned.  The  SYMAP1  steps  needed  include  a 
substitution  of  4  for  X,  a  substitution  of  -  3.25  for  Y,  and  a 
simplification.  The  substitutions  could  be  reversed  if  preferred,  and 
simplification  could  be  done  after  each  substitution  if  display  of 
intermediate  results  is  desired.  (See  Figure  8.) 

Substitution  of  constants  for  variables  in  trigonometric  and 
similar  expressions  is  possible,  but  simplification  is  not  complete  in 
such  cases  at  present.  Future  extensions  to  SYMAP1  may  provide  for 
additional  evaluations. 
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Given: 


X  100 

Y  101 

EXY  3.*X**.5*Y+X*Y**2. 


VX  4.$ 

VY  -3-25$ 


Steps : 


CONVERT 

PXY 

CP 

CSBEE 

VX 

X 

CSIME 

CPI 

CP2 

CSBEE 

VY 

Y 

CSIME 

CP3 

CPU 

RECONVERT 

CPU 

RP 

Results  in  SYMA.P1  notation: 


Results: 

P(X,Y) 

CP 

CPI 

P(U,Y) 

CP2 

CP3 

P(4,-3.25) 

cp  3.$ioo**.5^ioi**i4+i4ioo^i4,ioi**2.$ 

cpi  3  44  4** .  5$  ,ioi**i  4+i  4^  4**i  4,ioi**2  4 

CP2  64ioi**i4+44ioi**24 

cp3  64(-3.25$)**i4+^4(-3.25$)**24 

CPU  22.75$ 

RP  22.75 


Figure  8.  Evaluations 
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VIII.  VECTOR  OPERATIONS 


A  limited  number  of  operations  on  n-tuples  of  algebraic  forms  have 
been  provided  in  SYMAP1.  Others  could  be  added.  Currently  available 
are  VFRINT,  VCONVERT,  VKECONV,  VADD,  VSUBTR,  and  VSIME,  which 
respectively  print  an  n-tuple  of  consecutive  forms,  convert  an  n-tuple 
into  canonical  notation,  reconvert  from  canonical  form,  add  two  n-tuples, 
subtract  two  n-tuples,  and  simplify  two  n-tuples. 

In  each  case  n,  the  order  of  the  n-tuple,  is  specified  at  card 
column  4l  as  a  one-digit  positive  integer.  The  other  parameters  (at 
columns  11,  21,  31  as  needed)  are  those  which  would  be  needed  for  the 
first  form  of  the  n-tuple  in  non-vector  operations. 

Thus  VERINT  XI  b  b  7  would  cause  printing  of  the  7-tuple 
whose  elements  are  the  form  XI  and  its  six  successors,  whatever  their 
names.  The  various  forms  need  not  be  similar  with  respect  to  length  or 
any  other  feature .  Similarly  VCONVERT  ETA  IETA  b  3  would 

cause  ETA  to  be  converted  and  labeled  IETA  and  would  then  cause  the  next 
two  forms  after  ETA,  whatever  their  names,  to  be  converted  and  given 
names  much  like  IETA  by  adding  unity  for  each  in  the  9"th  character  of  the 
name.  VKECONV  acts  analogously,  as  does  VSIME. 

VADD  and  VSUBTR,  because  of  their  extra  parameter,  are  very 
slightly  different  from  the  above .  Thus  VADD  AB  CD  EF  4 

would  add  the  elements  of  the  4-tuple  whose  first  element  has  the  name 
AB  to  the  corresponding  elements  starting  at  name  CD  and  give  the 
(unsimplified)  results  the  names  EF  through  EFbbbbbb3 ,  respectively. 

These  vector  operations  merely  save  the  user  some  repetitive  card 
punching  when  very  similar  manipulations  are  needed  on  several  sets  of 
algebraic  forms.  The  user  must  be  careful  to  keep  his  sets  adjacent, 
as  it  is  the  adjacency  and  not  the  similarity  of  names  that  matters  in 
vector  operations.  Incidentally,  elements  of  n-tuples  may  be  manipulated 
as  independent  individuals  if  desired  by  using  their  individual  names 
in  the  usual  way.  (See  Figure  9*) 
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Figure  9-  Vector  Operations 


IX.  CONCHJSION 

The  present  SYMAP1  has  proven  to  be  useful  -within  its  realm  of 
application.  With  it  manipulations  can  be  carried  out  which  would  be 
extremely  tedious  if  done  manually,  and  the  likelihood  of  simple  errors 
of  arithmetic  is  reduced.  Various  alternatives  can  be  tried  with  little 
extra  effort.  Thus  a  researcher  can  be  freed  for  more  challenging 
investigations  and  more  quickly  be  shown  the  effects  of  changes  of 
attack  on  a  given  problem.  Potential  users  of  SYMAP1  are  encouraged  to 
try  it. 

Additional  features  will  be  added  to  the  program  as  they  are 
developed.  Several  of  these  are  planned  for  the  near  future,  and  still 
others  are  currently  being  studied.  Potential  users  of  SYMAP1  on 
BRLESC  1  are  encouraged  to  discuss  their  needs  with  the  author. 
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